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ABSTRACT 

We present a study of the evolution of the galaxy Velocity Dispersion Function (VDF) from z = to z = 1.5 
using photometric data from the UDS and NMBS COSMOS surveys. The VDF has been measured locally 
using direct kinematic measurements from the Sloan Digital Sky Survey, but direct studies of the VDF at 
high redshift are difficult as they require velocity dispersion measurements of many thousands of galaxies. 
Taylor et al. (2010) demonstrated that dynamical and stellar mass are linearly related when the structure of the 
galaxy is accounted for We show that the stellar mass, size and Sersic index can reliably predict the velocity 
dispersions of SDSS galaxies. We apply this relation to galaxies at high redshift and determine the evolution 
of the inferred VDF. We find that the VDF at z ^ 0.5 is very similar to the VDF at z = 0. At higher redshifts, 
we find that the number density of galaxies with dispersions < 200km/s is lower, but the number of high 
dispersion galaxies is constant or even higher. At fixed cumulative number density, the velocity dispersions 
of galaxies with log A^[ Mpc"-'] < -3 .5 increase with time by a factor of ^ L4 from z ^ L5 - 0, whereas the 
dispersions of galaxies with lower number density are approximately constant or decrease with time. The VDF 
appears to show less evolution than the stellar mass function, particularly at the lowest number densities. We 
note that these results are still somewhat uncertain and we suggest several avenues for further calibrating the 
inferred velocity dispersions. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: elliptical 
and lenticular, cD 



1. INTRODUCTION 

Stellar velocity dispersion is a fundamental property of 
galaxies. Through M ex R-^a\jG it provides a characteriza- 
tion of the galaxy mass and is a key axis in the Fundamental 
Plane of early-type galaxies (e.g., Djorgovski & Davis 1987; 
Bernardi et al. 2003). Further, it appears to correlate strongly 
with many other properties of galaxies, such as specific star 
formation rates, galaxy color, and black hole mass (Magor- 
rian et al. 1998; Franx et al. 2008; Trujillo, Ferreras, & de la 
Rosa 2011). 

The evolution of an individual galaxy's velocity dispersion 
carries information about the physical mechanisms responsi- 
ble for its growth. Broadly, processes that increase mass more 
efficiently than size, such as central gas accretion and the re- 
sulting star formation, could increase velocity while processes 
such as minor merging or mass loss in galactic winds could in- 
crease the overall size of a galaxy and decrease the velocity 
dispersion. The central dispersions may also be fairly stable 
with time and reflect the central dark matter potentials at the 
time when the galaxy started forming (Loeb & Peebles 2003). 
In this scenario, high redshift galaxies form the central re- 
gions of local massive galaxies and still retain signatures of 
their earliest progenitors. 

' Department of Astronomy. Yale University, New Haven, CT 06520- 
8101 

^ Sten'ewacht Leiden, Leiden University, NL-2300 RA Leiden, Nether- 
lands 

'European Southern Observatoiy, Alonso de CUrdova 3107, Casilla 
19001, Vitacura, Santiago, Chile 

* Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 

' Carnegie Observatories, Pasadena, CA 91101 

^ Hubble Fellow 

' Max-Planck-lnstitut fiir Astronomie, Konigstuhl 17, D-691 17 Heidel- 
berg, Germany 



Furthermore, central velocity dispersion lies at the inter- 
section of observational properties and quantities predictable 
from simulations. Observational evidence about the shape 
and evolution of the VDF can be directly compared to pre- 
dictions from cosmological simulations. Finally, masses of 
the supermassive black holes in the centers of galaxies dis- 
play the strongest correlations with host galaxy velocity dis- 
persions - understanding the distribution of velocity disper- 
sions would hold further significance for the study of black 
holes and AGN. 

In the local universe, the VDF has been measured directly 
(e.g., Sheth et al. 2003; Mitchell et al. 2005; Choi, Park, & 
Vogeley 2007). While measurements of stellar velocity dis- 
persion are possible out to z ^ 1, they become prohibitive for 
large samples at z > 1 . At z = we can predict velocity disper- 
sions of galaxies very well from their photometric properties. 
Here we present a new approach to describe the VDF based 
on photometric predictions and constrain its evolution out to 
z ^ 1.5. For this work we assume a concordance cosmology 
(//o = 70km/sMpc-i, f^M = 0.3 &rjA = 0.7). 

2. DATA 

2.1. Spectroscopic Redshifts, Velocity Dispersions, Sizes and 
Stellar Masses at Z'^Q 

The z ^ data is based on analysis of several publicly avail- 
able catalogs from SDSS DR7. Galaxies are included from 
0.05 < z < 0.07, using photometric information from the main 
DR7 catalogs (Abazajian et al. 2009) and redshift and ve- 
locity dispersions from the Princeton pipeline. The sample 
is selected to have good photometric measurements with the 
SDSS Science Primary flag, and low relative errors in velocity 
dispersion (< 10%). All velocity dispersions are aperture cor- 
rected to rj/8 using ctq = cra;,(8.0rap/re)"°^^ based on the best- 
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fit correction to the SAURON sample (Cappellari et al. 2006), 
where Vap = 1".5 is the radius of the SDSS spectroscopic fiber 
We adopt the best-fit Sersic (1968) effective radii in the r' 
band from the NYU-VAGC (Blanton et al. 2005). These sizes 
are based on fits to azimuthally averaged light profiles and 
are equivalent to circularized effective radii. Stellar masses 
are computed by the MPA-JHU group from the DR7 best-fit 
model magnitude^ (Brinchmannetal. 2004). We adopt these 
Mi^jL ratios and derive stellar masses from the luminosity of 
the best-fit Sersic model for all galaxies. 

2.2. Sizes and Stellar Masses at z> 

We use two samples of high redshift galaxies in this paper 
The first is the 0.77 deg^ UKIDSS Ultra-Deep Survey (UDS) 
K-selected galaxy catalog (Williams et al. 2009, 2010). This 
catalog includes near-infrared photometry (JHK) from the 
UKIDSS UDS (Lawrence et al. 2007; Warren et al. 2007), 
optical imaging from the SXDS (Sekiguchi & SXDS 2004) 
and 3.6/4.5/xm data from the SWIRE survey (Lonsdale et al. 
2003). Circularized effective radii (r^ = ^/ab) were measured 
in the J, H and K images for all bright sources (K < 22 A) 
from Sersic (1968) models fitted with GALFIT (Peng et al. 
2002). A full description of the catalog and size fitting as well 
as associated tests can be found in Williams et al. (2010). 

The second is the NEWFIRM Medium Band Survey 
(NMBS) of the COSMOS field, which we briefly summarize 
below. For an in-depth description of the survey see Whi taker 
et al. (2011). The 0.21 deg^ field includes medium-band 
NIR imaging (J1,J2,J3,H1,H2,K) from the May all 4m tele- 
scope, optical imaging (ugriz) from the CFHT Legacy Sur- 
vey and Subaru, IRAC imaging and MIPS 24 /im data. Sizes 
for the Cosmos galaxies are measured from the Ks WIRCam 
Deep Survey (WIRDS) 0".186/pixel images (Bielby et al., 
in prep), which have seeing of ^ 0".6-0".7 and the v. 1.3 
ACS F814W mosaic from Scoville et al. (2007). For de- 
tails of galaxy modeling see van Dokkum, et al., in prep. 
- we summarize as follows. Circularized effective radii, 
axis ratio and Sersic indices are determined for all galaxies 
brighter than K < 22 using Sersic models convolved with a 
position-dependent PSF with GALFIT (Peng et al. 2002). 
We test our procedure using independent high spatial reso- 
lution data from a small area of overlapping WFC3 imag- 
ing from the 3D-HST survey (van Dokkum et al., in prep). 
The interpolated WIRDS/ACS sizes are consistent with the 
WFC3 sizes down to ^ 0".15 with a biweight mean of 
log (re,ACS&wiRDs/'"f.wFC3) = 0".0009 and a scatter of 0.1 dex 
(Figure [T]). 

In both fields sizes are measured down to an arbitrary limit; 
however, as most of these sizes are based on ground-based 
imaging we assign a minimum size for all galaxies of r^ = 
0".2. 

Photometric redshifts for all galaxies in the UDS and COS- 
MOS surveys were calculated using the EAZY code (Bram- 
mer, van Dokkum, & Coppi 2008). Stellar masses were 
computed using FAST (Kriek et al. 2009), from Bruzual & 
Chariot (2003) stellar population synthesis models with So- 
lar metallicity and a Chabrier (2003) Initial Mass Function 
(IMF). 

3. INFERRED VELOCITY DISPERSION 



' http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 1 . — Sizes measured from WFC3 F140W imaging compared to sizes 
interpolated between measurements from ACS F8 1 4W and WIRDS Ks imag- 
ing. 

To derive a VDF at z > 0, we now describe how to estimate 
the stellar velocity dispersion of galaxies from their photomet- 
ric properties. According to the virial theorem, the dynamical 
mass of a galaxy is proportional to the product of the square 
of its measured dispersion and effective radius. 



Mdy„=K,alr,/G 



(1) 



Stellar mass has been shown to be proportional to dynam- 
ical mass for local galaxies (e.g., Taylor et al. 2010). How- 
ever, Taylor et al. (2010) demonstrated that simple estimates 
of dynamical mass based on homology exhibit residual trends 
with galaxy structural properties and introduced a more ro- 
bust structure-corrected dynamical mass, which we use to cal- 
culate inferred velocity dispersion based on photometric esti- 
mates of stellar mass and observed size. We adopt a Sersic de- 
pendent virial constant, Kv(n) (Bertin, Ciotti, & Del Principe 
2002) 



K„{n) = 



73.32 



10.465 + (n- 0.94)2 



; + 0.954. 



(2) 



In order to predict the central velocity dispersion of a galaxy 
based on its size and stellar mass, we define 



o-,„/ = 



GMdyn 



Ky(n)re V Ki,{n)re ' 



GM, 



(3) 



where K,{n) = K,{n)[-^^ 



Since Taylor et al. (2010) 



showed that Ky{n) ( -^ j depends only weakly on mass, we 

adopt the average ratio of stellar to total mass, calibrated such 
that the median amf = ctq for 2.0 < logcro < 2.4 in the SDSS. 

We find that (^)= 0.557, therefore KJn) = 0. 557 K,in). 

The central and inferred velocity dispersions agree well for 
the SDSS, with a 1 -u error of 0.06 dex (Figure|2^). This is re- 
markable, particularly when considering that no morphologi- 
cal selections are applied: after applying the Sersic dependent 
virial constant, this relation holds for both ellipticals and spi- 
ral galaxies. This implies that the velocity dispersion, ao, can 
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Fig. 2. — (a) Measured central velocity dispersion vs. inferred velocity dispersion of galaxies in the SDSS DR7. Galaxies are binned by logtr,-,,/- = 0.05 and 
the running mean and rms are shown in red for bins containing at least 100 galaxies. There is a very good correlation with a scatter of ~ 0.06 dex. High redshift 
(1 < z < 2.5) galaxies are included as large colored triangles, (b) Effective radii from NYU-VAGC best-fit Sersic models in the r' band vs. Stellar Mass in the 
SDSS, color coded by median measured velocity dispersion. Lines of constant inferred dispersion are included as dashed black lines. Both panels illustrate that 
the stellar mass and size of galaxies can be used to predict velocity dispersions. 



be well predicted from M^ and r^- Figure |2j5 shows the size- 
mass relation for SDSS, colored by median measured disper- 
sion for bins including a minimum of five galaxies. We use 
the median n{M^) to overplot lines of constant inferred dis- 
persion. Again, the measured and inferred dispersions agree 
quite well. 

Also included in Figure|2^ are published measurements for 
galaxies at 1 .0 < z < 2.5 (Cappellari et al. 2009; van Dokkum, 
Kriek, & Franx 2009; Newman et al. 2010; Onodera et al. 
2010; van de Sande et al. 2011) (large colored triangles), as- 
suming errors due to SED modeling of M,^ = 0.2 dex. The bi- 
weight mean offset in CTo/cr,„/ is small (-0.056 ±0.025 dex) 
and the high observed scatter of 0.131 ±0.028dex is fully 
consistent with scatter due to measurement errors. The re- 
lation does not appear to evolve significantly out to high red- 
shift, but we emphasize that current samples of high-z dynam- 
ical measurements are small and biased toward high disper- 
sions. Motivated by Figure |2^ we assume that the SDSS re- 
lations between velocity dispersion, stellar mass, effective ra- 
dius and Sersic index do not evolve with redshift and calculate 
inferred dispersions for the UDS and COSMOS catalogs. To 
limit bandpass-dependent effects we calculate cr,„y separately 
using UKIDSS J and K sizes for the UDS and ACS F814W 
and WIRDs Ks sizes for Cosmos. Next, we interpolate cr,,,/ 
to a rest-frame wavelength of 6200 , the central wavelength of 
the r' filter. We note that for apparently extremely compact 
galaxies with measured sizes < 0".2, these can be interpreted 
as minimum inferred velocity dispersions. 

4. VELOCITY DISPERSION FUNCTION 

A number of estimates of the z ^ VDF of early-type 
galaxies exist based on the SDSS (Sheth et al. 2003; Mitchell 
et al. 2005; Choi et al. 2007; Chae 2010). These studies are 
in reasonably good agreement and show that the z = VDF 
appears to be well fit by modified Schechter (1976) functions. 
Sheth et al. (2003) present a correction to the VDF for the 



contribution due to late-type galaxies, based on circular ve- 
locities. In the top row of Figure[3]we show both the fit to the 
early-type VDF (thin dotted black lines) and the total VDF 
(thick solid black lines), which includes the correction for 
late-type galaxies. This latter function is calculated by fitting 
a polynomial to the data in Figure 6 in Sheth et al. (2003). 

Measurement errors of velocity dispersions have a non- 
negligible effect on the shape of the inferred VDF, particularly 
on the steep, high dispersion tail. We model this effect by pro- 
ducing a mock catalog that intrinsically follows the total Sheth 
et al. (2003) VDF. In §3, we measured a scatter of 0.06 dex 
in inferred dispersion, which is comprised of intrinsic scatter 
in this relation combined with measurement errors in velocity 
dispersion. From the median relative error in velocity disper- 
sion for galaxies with g > 200 km/s (0.028 dex) we estimate 
that the intrinsic scatter in the relation between measured and 
inferred dispersion is 0.053 dex. The thick dashed black lines 
in Fig. [3] show the SDSS VDF when this scatter is added to 
the mock velocity dispersions. 

The top panel of Figure [3] shows the inferred VDF for each 
field in all redshift ranges and values are included in Table [T] 
Errors reflect Poisson errors in each velocity dispersion bin. 
Incompleteness is mostly caused by galaxies that are too faint 
to measure reliable sizes. We define the dispersion complete- 
ness limit (grey regions in Figure[3j for each field and redshift 
range as the dispersion at which the 95% completeness plus 
0.06 dex scatter about a linear fit to the K magnitude - inferred 
dispersion relation reaches the magnitude limit of the size cat- 
alogs. 

In the lowest redshift range, 0.3 < z < 0.6, the observed 
VDF matches the local VDF quite well in the region dom- 
inated by early-type galaxies, especially when compared to 
the scattered VDF. At low dispersions the z ^ 0.5 function 
is intermediate between the z = VDF for early-types only 
(black dashed line) and for all galaxies (black solid line). In 
the higher redshift bins, the VDF appears to flatten. Low 
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Fig. 3. — Top row. VDFs locally (black) as measured from SDSS and at higher redshifts for the two fields (color). Best-fit modified Schechter (1976) fits to 
the SDSS measured VDF are shown from Sheth et al. (2003) (thin dotted - early type VDF, thick solid - total VDF, thick dashed - total scattered VDF). Higher 
redshift number densities are based on inferred velocity dispersions assuming ^ = calibration. Individual colored lines reflect the measured VDF in each field 
(solid for UDS and dashed for Cosmos) and the circular data points show the volume-weighted average. Bottom row: Cumulative VDFs for each redshift bin. 
Colors and symbols are the same as above. Three cuts in number density {horizontal gray lines) highlight three possible evolutionary behaviors. Galaxies at 
the lowest number density threshold, log ^[Mpc"^] < -5.5, increase in velocity dispersion from a redshift z ~ 1.5 to z ^ 0. Galaxies at a number density of 
log ^[Mpc"^] < —4.5 appear to be near the crossover and experience Httle evolution in velocity dispersion to z ~ 0. Finally galaxies at the highest number density 
threshold of log ^[Mpc"^] < -3.5 appear to decrease in velocity dispersion with time. 



dispersion galaxies become less common and high disper- 
sion galaxies more common, with an approximately constant 
crossover point. The colored points show results for the UDS 
(diamonds) and COSMOS (triangles) fields separately. The 
flattening of the VDF with redshift is seen in both fields, sug- 
gesting that it is not due to field-to-field variations. We note 
that the COSMOS field has a cluster at z=0.7 and a overden- 
sity at z ~ 1.0 (Whitaker et al. 2011). These structures are 
reflected in the higher normalization of the COSMOS VDF in 
the second and third redshift bins. 

5. CUMULATIVE VELOCITY DISPERSION FUNCTION 

We also consider the evolution of velocity dispersions of 
galaxies at constant cumulative number density. Assuming 
that the rank order of galaxies remains the same, a num- 
ber density selection selects the same population of galax- 
ies as a function of redshift (e.g.. Wake et al. 2008; van 
Dokkum et al. 2010). We show the cumulative VDF in the 
bottom row of Figure [3] and select galaxies at three num- 
ber densities: log A^[Mpc"^^] = -3.5,-4.5,&-5.5 (horizontal 
gray lines). The crossing point between each number density 
threshold and the cumulative VDF evolves differently with 
redshift. 

In Figure IH we show the evolution of the velocity disper- 
sion at fixed cumulative number density {(Jcross) with redshift. 
The rare high-dispersion galaxies (log A^[Mpc"^] < -5.5) ex- 
hibit a weak trend with even higher dispersions at high red- 
shift and a slope of 0.066 ± 0.059. While it is difficult to 
decrease a galaxy's velocity dispersion with time, it has been 
suggested that this can occur via extensive minor merging. 



which can efficiently increase size while slowly increasing 
mass (e.g., Bezanson et al. 2009; Hopkins et al. 2009; van 
Dokkum et al. 2010). An abundance of evidence has demon- 
strated that at a given mass, galaxies have smaller sizes at high 
redshift than in the local universe (e.g., Daddi et al. 2005; 
Trujillo et al. 2006; van Dokkum et al. 2008; Cimatti et al. 
2008; van der Wei et al. 2008; Franx et al. 2008; Damjanov 
et al. 2009). We emphasize that these are the most extreme 
galaxies in the universe - with inferred velocity dispersions of 
> 400 km/s at z ^ 0. The fact that this trend is weak sug- 
gests that the structure correction to the inferred dispersions 
may mitigate the evolution in density within r^. At the higher 
threshold of log A^[Mpc"^] < -4.5 galaxies appear to have al- 
ready attained the central velocity dispersions that they will 
have at z = 0, with a slope of 0.014 ± 0.054. If these galax- 
ies maintain constant central dispersion while continuing to 
evolve from z ^ 1 .5 to the present it could be evidence that 
they are undergoing inside-out growth in which the central 
dynamics of the stars are set by some early processes and the 
overall structure of the galaxy builds up around the core or 
that they undergo size growth at the same rate as mass growth. 
Finally, for the galaxies above the highest density threshold of 
logA^[Mpc"^] <-3.5 velocity dispersion increases with time. 
In this case the negative slope of -0.098 ± 0.053 is signifi- 
cant. This may suggest the importance of gas accretion and 
star formation, which could increase the stellar masses of the 
galaxies. 

We note that the dispersion function appears to evolve 
slower than the stellar mass function: stellar mass increases 
by a factor of ^ 2 for the most extreme galaxies (log Mi, > 11) 
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TABLE 1 
Inferred and Cumulative Velocity Dispersion Functions 



logo-,,,/ [km/s] 



Redshifts 



2.0-2.1 



2.1-2.2 



2.2-2.3 



2.3-2.4 



2.4-2.5 



Field 



log<J>[Mpc ^dex 



2.5-2.6 



2.6-2.7 



0.3 <; 


:<0.6 


UDS 
Cosmos 


-2.25 ±0.02 
-2.24 ±0.04 


-2.27 ±0.02 
-2.27 ±0.04 


-2.38 ±0.03 
-2.28 ±0.04 


-2.62 ±0.04 
-2.32 ±0.05 


-3.32 ±0.08 
-2.83 ±0.08 


-4.47 ±0.31 
-3.69 ±0.22 


<-4.77 
<-4.29 






Total 


-2.24 ±0.02 


-2.27 ±0.02 


-2.36 ±0.02 


-2.52 ±0.03 


-3. 14 ±0.06 


-4.11 ±0.18 


< -4.89 


0.6 <; 


:<0.9 


UDS 
Cosmos 


-2.44 ± 0.02 
-2.38 ±0.03 


-2.51 ±0.02 
-2.36 ±0.03 


-2.54 ±0.02 
-2.26 ±0.03 


-2.77 ±0.03 
-2.30 ±0.03 


-3.30 ±0.06 
-2.48 ±0.04 


-4.76 ±0.31 
-3.26 ±0.09 


< -5.06 
-4.58 ±0.43 






Total 


-2.43 ±0.02 


-2.46 ±0.02 


-2.45 ±0.02 


-2.60 ±0.02 


-2.92 ±0.03 


-3.82 ±0.09 


-5. 19 ±0.43 


0.9 <; 


:< 1.2 


UDS 
Cosmos 


-2.62 ±0.02 
-2.54 ±0.03 


-2.70 ±0.02 
-2.62 ±0.04 


-2.78 ±0.03 
-2.63 ±0.04 


-2.94 ±0.03 
-2.69 ±0.04 


-3.40 ±0.05 
-2.93 ±0.06 


-4.06 ±0.12 
-3.45 ±0.10 


-4.91 ±0.31 
-3.82±0.15 






Total 


-2.60 ±0.02 


-2.68 ±0.02 


-2.74 ±0.02 


-2.86 ±0.03 


-3.23 ±0.04 


-3.81 ±0.08 


-4.33 ±0.14 


1.2 <; 


:< 1.5 


UDS 
Cosmos 


-2.65 ±0.02 
-2.94 ±0.05 


-2.64 ±0.02 
-2.92 ±0.05 


-2.71 ±0.02 
-3.11 ±0.06 


-2.87 ±0.03 
-3.11 ±0.06 


-3. 17 ±0.04 
-3.53±0.10 


-3.91 ±0.09 
-3.76±0.13 


-4.68 ±0.22 
-4.33 ±0.25 






Total 


-2.70 ±0.02 


-2.69 ±0.02 


-2.78 ±0.02 


-2.92 ±0.02 


-3.23 ±0.04 


-3.87 ±0.07 


-4.57 ±0.16 



logo-,,,/ [ki-n/s] 



Redshifts 



>2.0 



>2.1 



>2.2 



>2.3 



>2.4 



>2.5 



Field 



logA'[Mpc 



>2.6 



0.3<z<0.6 


UDS 
Cosmos 


-1.74 ±0.01 
-1.64 ±0.02 


-1.90 ±0.02 
-1.77 ±0.02 


-2.15 ±0.02 
-1.93 ±0.03 


-2.53 ±0.03 
-2. 19 ±0.04 


-3.29 ±0.08 
-2.77 ±0.08 


-4.47 ±0.31 
-3.69 ±0.22 


<-4.77 
<-4.29 




Total 


-1.71 ±0.01 


-1.87 ±0.01 


-2.09 ±0.02 


-2.42 ±0.03 


-3.09 ±0.05 


-4.11±0.18 


< -4.89 


0.6<z<0.9 


UDS 
Cosmos 


-1.93 ±0.01 
-1.64 ±0.01 


-2.08 ±0.01 
-1.73 ±0.02 


-2.29 ±0.02 
-1.84 ±0.02 


-2.65 ±0.03 
-2.05 ±0.02 


-3.28 ±0.06 
-2.41 ±0.04 


-4.76 ±0.31 
-3.24 ±0.09 


< -5.06 
-4.58 ±0.43 




Total 


-1.84 ±0.01 


-1.96±0.01 


-2.13 ±0.01 


-2.41 ±0.02 


-2.87 ±0.03 


-3.81 ±0.09 


-5. 19 ±0.43 


0.9 <z< 1.2 


UDS 
Cosmos 


-2.11 ±0.01 
-1.94 ±0.02 


-2.28 ±0.01 
-2.07 ±0.02 


-2.48 ± 0.02 
-2.22 ±0.02 


-2.78 ±0.03 
-2.43 ±0.03 


-3.30 ±0.05 
-2.78 ±0.05 


-4.00 ±0.11 
-3.30 ±0.08 


-4.91 ±0.31 
-3.82±0.15 




Total 


-2.06 ±0.01 


-2.22 ±0.01 


-2.40 ±0.01 


-2.66 ±0.02 


-3. 10 ±0.03 


-3.70 ±0.07 


-4.33±0.14 


1.2 <z< 1.5 


UDS 
Cosmos 


-2.06 ±0.01 
-2.35 ±0.03 


-2. 19 ±0.01 
-2.48 ±0.03 


-2.38 ±0.02 
-2.68 ±0.04 


-2.66 ±0.02 
-2.89 ±0.05 


-3.08 ±0.03 
-3.29 ±0.08 


-3.84 ±0.08 
-3.66 ±0.12 


-4.68 ±0.22 
-4.33 ±0.25 




Total 


-2.12 ±0.01 


-2.25 ±0.01 


-2.44 ±0.01 


-2.71 ±0.02 


-3.13±0.03 


-3.79 ±0.07 


-4.57±0.16 



from z ^ 2-0, with even stronger evolution at lower masses 
(Brammer et al. 201 1). 

6. DISCUSSION 

We stress that all of these results are preliminary and rest on 
a multitude of assumptions. First of all, our samples of galax- 
ies with size measurements may suffer from incompleteness, 
especially at the low mass limit. We have tried to place thresh- 
olds above which our sample should be complete, however 
these limits may be underestimated. 

The majority of galaxy sizes are measured from ground- 
based imaging, which could be uncertain, especially for com- 
pact galaxies at high redshift. In Figure[T]we have shown that 
the WIRDS/ACS sizes agree quite well with those measured 
from WFC3, exhibiting a scatter of only 0.1 dex. However, at 
the high dispersion end of the VDF, the accuracy of individ- 
ual size measurements becomes more important than overall 
trends. Future studies would benefit from careful fitting of 
profiles exclusively based on space-based imaging. This will 
also improve the accuracy of our Sersic indices. 

The most tentative assumption remains that the relation be- 
tween inferred and measured velocity dispersions does not 
evolve. There are many factors that contribute to this relation. 



Aperture corrections depend on galaxy profiles and therefore 
any systematic evolution in the light distribution of galaxies 
would affect the measured velocity dispersions. However, this 
is a relatively weak effect as the aperture corrections are small. 
More important may be the dependence ofKy on Sersic index, 
which could evolve if galaxy structure evolves in complex 
ways. Furthermore, our calibration of ^*(n) depends strongly 
on the average ratio of stellar to total dynamical mass, which 
may evolve with time and depend on the stellar mass of the 
galaxy. Taylor et al. (2010) found that M^/Mdy„,„ still ex- 
hibits trends with stellar mass at fixed «, suggesting that even 
at z = the models for Ky{n) do not completely describe the 
non-homology of galaxies. Any of these effects could pref- 
erentially influence certain types of galaxies, introducing sys- 
tematic trends in the shape of the ctq vs. (Ti„f relation. The 
scatter in the relation may also evolve with time; this could 
explain the apparent increase in number of galaxies with high 
inferred dispersions at z > 0.6. Finally, we assume a constant 
Chabrier (2003) IMF, which could also differ with redshift or 
halo mass (e.g., van Dokkum & Conroy 2010). 

Measuring the entire VDF using dynamical measurements 
would require an extremely large spectroscopic survey. A 
more economical approach is to better calibrate the evolution 
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Fig. 4. — Velocity dispersion at tlie three number density tliresliolds 
(log A'[ Mpc"^^ ] = -3.5, -4.5, &- 5.5) shown in the bottom row of Figure [3] 
Open symbols represent the unscattered (small symbols) and scattered (large 
symbols) Sheth et al. (2003) total VDF Filled symbols represent rjcmss for 
the higher redshift VDFs. Vertical error bars include the 0.06 dex error in in- 
ferred velocity dispersion divided by the square root of the number of galaxies 
interpolated between two nearest bins in logcr from Figure[3] We include an 
additional error of 0.05 dex (error bar in lower left comer), added in quadra- 
ture to all points, to encompass additional sources of error. 

of the inferred VDF. With a large sample of uniformly se- 
lected galaxies spanning a large range of velocity dispersions 



and lower measurement errors, one could characterize a-mf as 
a function of stellar mass, Sersic fits and redshift. 

Finally, we note that Chae (2010) measured the VDF in a 
similar redshift range (z ^ 0.3 - 1 .0) from 30 lensing galaxies 
and found quite different results. Chae finds that at the low 
dispersion end (logcr < 2.3) the VDF of galaxies agrees quite 
well with the local VDF, and that it steepens at higher red- 
shift. At the high dispersion end, the Chae (2010) results are 
based on a lack of arcs with large separations. It is possible 
that this sample of lensing galaxies preferentially misses such 
lenses. Alternatively, it is possible that the intrinsic VDF does 
fall below the local VDF at high dispersions, but the scatter 
increases rapidly with redshift, producing the opposite evolu- 
tion in our study. Further calibration of the inferred velocity 
dispersions at higher redshift is necessary to probe whether 
the two results are inconsistent. 

Although measurements of the VDF at higher redshift may 
never be as precise as those measured from the SDSS, study of 
the evolution of the inferred velocity dispersion function using 
more accurate and redshift-dependent dynamical calibrations 
will be a powerful tool for studying the growth and evolution 
of galaxies and the halos that they occupy. 

We thank the referee for his/her feedback and acknowledge 
support from NSF grant AST-0807974 and HST grant GO- 

12177.01A. 
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